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1. Introduction 

Currently developed mean-field theories based on energy density functionals ambitiously 
aim to properly explain a large range of nuclear properties and in many fields they 
have been quite successful. In this paper I present the results of applying of the 
UNEDFO functional [lj to describe low-energy collective excitations in the chain of 
even-even 118 “ 144 Xe isotopes. The treatment of collective properties is based on the 
Adiabatic Time Dependent HFB (ATDHFB) theory, which leads to a construction of 
a collective Hamiltonian from a microscopic, mean-field, input. More details of the 
applied methods can be found in [2][3|. An alternative approach to collective phenomena 
within a microscopic theory with phenomenological interactions employs the Generator 
Coordinate Method (GCM), see e.g. jU-Ej. One should also mention an attempt to 
describe collective phenomena, in the RPA context however, starting form realistic 
interactions and using the Unitary Correlation Operator method (UCOM) [TIE] • 

For many years Xe nuclei have been subject of extensive experimental and 
theoretical studies. Many of them were focused on a phenomenon of the double /3 
decay (confirmed experimentally in the case of 136 Xe) but there is also a rich literature 
concerning the collective properties related to changes of nuclear deformation. The long 
chain of Xe isotopes offers a good opportunity to study the evolution of these properties 
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as the number of neutrons increases as well as the role of nonaxiality. Let me mention 
only a few works which studied one or more of the even-even Xe isotopes: papers 
employing geometrical concepts and the Bohr Hamiltonian |2ll9LH2] . papers based on 
the interacting boson model [T3lfT6] . papers using truncated shell model space HM9]. 

Section [2] presents some basic facts on the general Bohr Hamiltonian, on methods 
which allow for its derivation from the mean-field theory and on the UNEDFO density 
energy functional. In section [3] I show the results of calculations concerning several low- 
spin energy levels, some E2 transitions in the 118 ~ 144 Xe nuclei as well as the comparison 
with experimental data. 

2. Theory 

2.1. Quadrupole variables, the Bohr Hamiltonian 

A consistent description of nuclear vibrational and rotational excitations as well as of 
possible couplings between them requires the use of quadrupole collective variables i.e. 
of the second rank with respect to the SO(3) rotation group. Such variables can be 
chosen in various ways, e.g. as parameters describing the shape of a nucleus [ 20 ll 2 Tj 
or the shape of a phenomenological one-particle potential (gum Eg. Within a self- 
consistent mean field theory such quadrupole variables 07 (in the laboratory frame) are 
chosen so as to be proportional to the components of the quadrupole mass tensor 

A 

02 , ~ (i) 

2—1 

where $ is a microscopic nuclear wave function which can be obtained by using effective 
interactions of the Skyrme [2] or the Gogny type [231124] or in the relativistic framework 
(RMF) [2511263 - The quadrupole variables can be equivalently expressed in the intrinsic 
frame (also called principal axes frame) by two deformation variables /3, 7 and three 
Euler angles (fl) describing the relative orientation of the laboratory and intrinsic frame. 
The j3, 7 variables are given by mean values of the operators Qo = Yhi=i^ z< i ~ r i ) an d 
Q 2 = 1 >/3 (xf - y'f) as follows 

/3 c osy = cq 0 , qo = ($|Q 0 |$) (2) 

/5 sin 7 = cq 2 , q 2 = ($|<3 2 |$) (3) 

with a conventional factor c = \ Jn/5/Ar 2 where r 2 = 3/5vqA 2 ^ 3 , vq = 1.2 fm. 

One should keep in mind that in some theoretical approaches, e.g. in the geometrical 
collective (Frankfurt) model [[27] . the deformation variables /3, 7 do not have a direct 
relation to a nuclear shape or mass distribution. Within a framework of the interacting 
boson model [28j the /3, 7 variables introduced by means of the so called coherent states 
are related rather with properties of valence nucleons and not of a spatial distribution 
of a nuclear density. 

The general properties of the quadrupole collective space as well as of functions and 
operators depending on the quadrupole variables can be found e.g. in [3]. The most 
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important, from the point of view of physical applications, is a Hamiltonian which we 
call the General Bohr Hamiltonian (GBH) and which can be expressed in the intrinsic 
frame as 

-f^Bohr — -hyib + T rot + V (4) 


T v ib = 


1 r 


wr \ /3 4 


dr- 


-B^)d p -d^(/3\l-Bg^) d. 


M I 


+ 


/3 sin37 

1 


a 


7IV ... sill3 7 5 /3 7 ) dp + "jfd'y ( \j — mn&yBpp ) d. 


T rot = - Y, 4 2 (^)/TG Jk = 4 B k (/3, 7 )/3 2 sin 2 ( 7 - 2nk/3) 


k =1 


where w = BggB^ - B 2 ^, r = B x B y B z 


(5) 

( 6 ) 
(7) 


The operators Ik(Tl),k = 1,2,3 are components of the angular momentum in the 
intrinsic frame. The Hamiltonian (j4j) contains seven functions that depend on 
deformation variables: the potential energy V and six functions Bgg , Bg- p H 77 , Bk, called 
mass parameters or inertial functions. 

One possible way to determine these seven functions consists in assuming for them a 
’reasonable’ form with some free parameters which are determined through comparison 
of calculated and experimental collective properties. 1 use another approach which is 
based on the ATDHFB (Adiabatic Time Dependent HFB) theory and which aims at 
calculating these functions starting from a microscopic theory. In this approach one does 
not introduce any additional free parameters and the prediction of collective properties 
is based solely on the knowledge of effective nucleon-nucleon interactions. 


2.2. The ATDHFB mass parameters 

In the following discussion it is assumed that the time evolution of a system is 
determined through a time dependence of several collective variables qg (not necessarily 
the quadrupole ones from eqs. (EHSJ)). The ATDHFB theory, based on an assumption of 
low collective velocities, gives an expression 4 j Bkjfkqj which is bilinear in velocities 
and which defines a metric tensor in the collective space. In the next step this 
expression is used to calculate the Laplace-Beltrami operator which is taken (up to 
the H 2 factor) as a kinetic energy part of a collective Hamiltonian. Functions Bkj (mass 
parameters) depend on collective variables. More details on the ATDHFB theory and 
mass parameters can be found in e.g. |3j[29] and papers cited therein. Below I briefly 
sketch some steps and give some formulas which are needed to calculate the general 
Bohr Hamiltonian starting from the UNEDFO energy functional. 

The so called cranking approximation ignores the Thouless-Valatin terms so that 
the mass parameters can be conveniently expressed through derivatives of a generalized 
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density matrix R(q) corresponding to the HFB state 4>(q). The derivative d qk R in the 
quasi-particle basis (in the doubled space) has the form 


■£~r) =F k = 

Qk J quasipart 



and the mass parameters read 

R Ff fj^ufk, 


fils 




fJ,U 


E a + E v 


f k = -n 


( 8 ) 

( 9 ) 


where E^ v are quasi-particle energies. If the matrix d qk R is known in a fixed single- 
particle basis the matrix F k can be calculated as 


F„ = B(d qh R ) 


fixed, doubled 


B 


+ 


( 10 ) 


where B is the Bogolyubov matrix for R 


B 


U+ V + \ 

v T u T J ■ 


( 11 ) 


Sometimes the following alternative expression for f k is useful 


fk,nu = ($| a v ay. \d qk §) (12) 

where are quasi-particle annihilation operators. 

In the case of quadrupole variables one obtains the deformation-dependent HFB 
state by constrained HFB calculations 


S ( l I l |// l(u , :1 1'I') = 0 with ('!' Qjir = q r l = 112 


(13) 


Then, it is easier first to discuss the vibrational mass parameters B qiqk ,k = 0,2 
from which Bpp, -Bg 7 , B^ n in the formulas ([5HTJ) can be calculated by a simple change 
of variables. The required derivatives d qk R should be calculated by numerical 
differentiation (see 01291 ED]) but most often one resorts to another (so called 
perturbative) approximation which relates the derivatives of the generalized density 
matrix to derivatives of the induced one-body Hamiltonian [3j. The constraints (eq fldjl 
lead to the extra term —A jQj in the induced one-body Hamiltonian and one can easily 
calculate a derivative with respect to Ay 




( w i) ^ 


E^ + EX 


(14) 


where 


Wj = U + {Qj) fixed- (U + (Qj) &xed V*) T 


Then, the derivatives d qk R are calculated using the relation 


_d_ 

dq k 


j 


d\j d 
dq k d\j 


(15) 


( 16 ) 
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and finally the derivatives d\j/dq k are obtained by inverting the matrix dq k /d\j, which 
can also be expressed through Wj, eq (1T5T) 

n _ k)j) 




Re E 


fw 


En + E u 


The moments of inertia are given by the Inglis-Belyaev formula 


J k = Y 


I [u + (j. 


k ) fixed 


v* 


(U+( Jk ) &xed V*) T ] 


fLD I 


E u + E u 


(17) 


(18) 


where (jfc)fixed is a matrix of the microscopic total angular momentum. 

In the case of the BCS approximation when the canonical basis is used formulas 
(TT511TB|) can be transformed to a simpler and better known form but in the general HFB 
approach (fT5l IT8|) are more useful. 


2.3. The UNEDFO Energy functional 

To construct the mean-field configurations I used the UNEDFO energy density 
functional, which is one of the results of a large scale project UNEDF |31j . The 
functional is described in detail in [I] and here I will elaborate only on some of its 
distinctive features. In the particle-hole channel UNEDFO is a ’standard’ Skyrme-type 
functional [32] with the spin-orbit term treated as in the Ski parametrization [33]. The 
pairing interaction is modelled as a sum of the standard (volume) plus density-dependent 
surface peaked 5 interaction 

<J(ri-r 2 ), q = n,p; p 0 = 0.16 frn -3 (19) 

and the Lipkin-Nogami method is used to avoid the pairing ’collapse’ for magic nuclei 
and their neigbours. The pairing strengths for protons and neutrons, Vff ,v are fitted 
simultaneously with other parameters determining the functional. A truncation of the 
quasi-particle space, required due to a zero-range of the pairing interaction is fixed by 
the condition for quasi-particle energies E jt < 60 MeV. Because it is well known that 
the ATDHFB mass parameters are quite sensitive to a diffuseness of the occupation 
number distribution I shall now present more details on the treatment of the pairing 
part of UNEDFO. 

The binding energies of the considered Xe isotopes are reproduced quite well by the 
UNEDFO functional. The RMSD (root mean square deviation) for 16 nuclei is equal 
to 0.454 MeV with the largest error B t h — B exp = —0.71 MeV for the 134 Xe isotope. 
The chain of isotopes contains 136 Xe with a magic number N = 82 of neutrons but it 
appears that due to the Lipkin-Nogami (LN) prescription the changes of the pairing 
properties along the chain are quite smooth. This can be seen in figure [1] where I plot 
the neutron and proton pairing energy vs the mass number. In addition I show a plot 
of the quantity A av + A 2 ,ln, where A av = ^ u^v^A^/ ^ and A 2 ,ln is a coefficient 
determined in the LN method. This quantity can be treated as an estimation of the 
pairing gap within the LN method, for more details see [ 34] , 


V q 

v o 


1 + 1 - 


P( r ) 
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Figure 1 . Plot of the pairing energy (upper panel) and of the effective pairing 
gap A av + \ 2 ,ln calculated for protons (♦)and neutrons (O ) at the deformation 
corresponding to a minimum of the potential energy. 


In conclusion I want to mention two newer functionals UNEDF1 [35j and UNEDF2 
132,135] which were constructed by extending the empirical dataset used in the fitting 
procedure. In the case of the UNEDF1 functional new data on a few fission isomers was 
added while for the UNEDF2 several single-particle level splittings were additionally 
considered. However, the RMSD for binding energies is significantly lower (around 
1.4 MeV) for UNEDFO than for UNEDF1 and UNEDF2 (around 1.9 MeV) hence the 
UNEDFO functional seems to be a good choice for a pilot study of collective properties 
in the region of medium-heavy nuclei. A further detailed study on the consequences of 
UNEDF1 and UNEDF2 for collective nuclear properties is currently in progress. 

3. Results of calculations, comparison with experiment 

The values of inertial functions and potential energy which enter the General Bohr 
Hamiltonian were calculated at 144 points forming a regular grid in the sextant 
(0 < f3 < 0.65) x (0 < 7 < 60°) in the deformation plane. The distance between 
the points is 0.05 and 6 ° in the fd and 7 directions, respectively. The mean-held wave 
functions were obtained using the code HFODD ver. 2.49t, see ['3Tj and references 
therein. 

3.1. Potential energy surfaces 

As can be seen in figure [ 2 ] there are three nuclei 134-138 Xe with a spherical minimum of 
the potential energy. Others exhibit deformed minima with (3 m i n in the range 0.1 — 0.25, 
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mostly on the prolate axis except for 118 Xe and 128 Xe which have slightly nonaxial 
minima with 7 min = 8 ° and 12°, respectively. The depths of the minima (relative to a 
spherical shape) are less than 2 MeV. In figures [3] and [4] I show full plots of the potential 
energy on the deformation space for a representative sample of four isotopes. One can 
notice a rather weak dependence of the potential energy on the 7 variable (7 softness), 
especially for lighter isotopes. 



A 



MeV 

-0.5 

“I- 

-1-1-1- 

—i—i—^ ♦ T T w. —i—i—r 

♦ ♦ 

♦ 

(b) 

-1.0 

♦ 


f V -V u ♦- 

v min v spher 


-1.5 

"i 




-2.0 

J_ 

_1_1_1_ 

_1_1_1_1_1_1_1_1_1_L 


114 116 118 120 122 124 126 128 130 132 134 136 138 140 142 144 


A 



A 


Figure 2. a) Value of the deformation /3 m ; n (left panel) and 7 m i n (right panel) at the 
minimum of the potential energy, (b) Depth (relative to a value at a spherical shape) 
of minima of the potential energy for the 114-144 Xe nuclei. 
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Figure 3. Plot of the potential energy (relative to a spherical shape value) for the 
118 Xe and 126 Xe isotopes. 
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Figure 4. Plot of the potential energy (relative to a spherical shape value) for the 
134 Xe and 144 Xe isotopes. 


3.2. Collective energy levels 

Having calculated the potential energy and mass parameters I performed a numerical 
diagonalization of the resulted Bohr Hamiltonian using the method described in [2, IDJ - 
The obtained eigenvalues can be directly compared with excited energy levels of positive 
parity and the corresponding collective wave functions can be then used to calculate 
matrix elements of various operators, in particular of the operators of electromagnetic 
E2 transitions. 

1 should add that all mass parameters (vibrational and rotational) were multiplied 
before the diagonalization of the Bohr Hamiltonian by a constant factor 1.3. The 
commonly quoted reasons for introducing such a factor refer to a simulation of the 
effects of including the Thouless-Valatin terms in the ATDHFB method or/and the 
effects of the so called pairing vibrations, see e.g. [2|l3] l23] l38 ] l39]. Some rather crude 
estimations of these effects give the value of the factor in the range 1.2 — 1.5. However, 
due to a lack of sufficiently quantitative calculations this factor must be treated as an 
additional parameter of the theory. Before presenting the results for the whole chain of 
Xe isotopes 1 will show consequences of introducing the scaling factor for energy spectra 
and B(E2) probabilities in the case of 126 Xe. Figure [5] contains plots of bands built on 
Oi, 2 2 and 0 2 levels. There are two sets of theoretical results: obtained with the scaling 
factor (sc = 1.3) and without the scaling (sc = 1). One can see that the scaling produces 
a ’shrinking’ of the spectra leaving a general picture similar in both cases. In addition, 
one can see that the scaling leads to a better agreement with experimental data (showed 
in figure [5] as well). A sample of B(E 2) results (theoretical with and without scaling 
and experimental) is shown in figure [6] The sample contains cases with both good and 
worse agreement between theory and experiment. It can be seen that the effect of the 
scaling the mass parameters is much smaller on the B(E2) probabilities than on the 
values of level energies. 

Then I compare theoretical energies of several low lying low spin levels 
(2 1? 4 1; 2 2 , 0 2 , 0 3 ) with experimental data [TO] for the considered chain of Xe isotopes. 
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Figure 5. Energy levels in 126 Xe. Comparison of theoretical results obtained using the 
scaling of mass parameters (sc = 1.3), without such scaling (sc = 1) and experimental 
data. 


126 Xe, B(E2) [W.u.] 



Figure 6. Selected E2 transitions in the 126 Xe nucleus. Theoretical results with the 
scaling factor (sc = 1.3), without one (sc = 1) and experimental data |16l are shown 
with (D),(A) and (• ), respectively 


These levels were chosen because of their role in analysing the band structure of nuclear 
spectra. 

One can conclude from plots in figures ITlfTTl that in general the theoretical results 
are in good agreement with experimental data, especially for the lighter part of the 
isotope chain (up to A = 130). One should also keep in mind that I do not fit any 
parameters to collective properties. Some significant discrepancies in the vicinity of 
N = 82 number of neutrons are not unexpected because the ATDHFB theory tends to 
perform better for more collective nuclei, i.e. with a larger number of valence nucleons 
(let us recall that there are only four valence protons in Xe isotopes). This effect is 
connected with an assumption of adiabatic motion of all nucleons in the varying mean- 
field. This assumption is strongly affected by a presence of closed shells. 
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A 

Figure 7. Experimental (■) [50] and theoretical (□) energy of the 2f level in the 
H4 _ i4 4 x e nuclei. 



A 

Figure 8. Same as in figure [TJ but for the 4f level. 



A 


Figure 9. Same as in figure 0 but for the 2f level. 
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A 

Figure 10. Same as in figure [7] but for the 0^" level. 



A 

Figure 11. Same as in figure [3 but for the O 3 " level. 


3.3. E2 transitions 

A detailed analysis of experimental data on electromagnetic transitions can provide 
important information about excited levels, see e.g. [U]. In the case of quadrupole 
excitations the most important are E2 transitions which are described by the collective 
operator 

Qtr l ’\P, 7) = (20) 

1 present the results of calculations of the B(E2) reduced transition probabilities for 
transitions 2\ —> Oi and 4i —> 2i as well as their comparison with evaluated experimental 
data from [40] in figures [13 and [£2J Again one can see that theoretical calculations 
reproduce the general behaviour of B(E2) quite well even as no free parameters (e.g. 
effective charges) were used. 

In the case of some Xe isotopes there are much more extensive experimental data 
on E2 transitions, see e.g. for 126 Xe [16] and for 128 Xe [421143] but I postpone a discussion 
of them to a subsequent publication. 
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A 

Figure 12. Experimental (• ) [ 30 ] and theoretical (□) B(E2) probability (in W.u.) 
for the transition 2f —»• Of level in the 114-144 Xe nuclei. 



A 

Figure 13. Experimental (• ) [40] and theoretical (□) B(E2) probability (in W.u.) 
for the transition 4f —> 2f level in the 114_144 Xe nuclei. 


4. Conclusions 

The paper presents the results of the first attempt to apply the UNEDFO energy 
functional to the theory of a nuclear collective motion. A correct description of the 
properties of the long chain of Xe isotopes considered is a demanding challenge to 
the theory, in particular to a framework with no free parameters that could be fitted 
to experimental data on collective levels. It can be argued that the results shown in 
section [3] are quite satisfactory and reproduce well the general tendencies seen in the 
energy spectra as well as E2 transitions in the Xe isotopes, with some exceptions around 
the semi-magic 136 Xe isotope. This contribution contains only a part of the obtained 
theoretical results, a more detailed analysis of the 120-128 Xe nuclei is in preparation. 
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